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Abstract 



We analyze doping dependent spectral intensities and Fermi surface maps 
obtained recently in Nd2- x Ce a; Cu04± ( 5 (NCCO) via high resolution ARPES 
measurements, and show that the behavior of this electron-doped compound 
can be understood as the closing of a Mott (pseudo) gap, leading to a quan- 
tum critical point just above optimal doping. The doping dependence of 
the effective Hubbard U adduced by comparing theoretical and experimental 
spectra is in resonable accord with various estimates and a simple screening 
calculation. 
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Angle-resolved photoemission spectra (ARPES) follow a remarkably different route with 
doping in the hole-doped La 2 _ :r Sr, r Cu04 (LSCO) |] compared to the electron- doped NCCO 
0. Starting from a Mott insulator near x = 0, LSCO exhibits the appearance of dynamic 
or static stripes [p]], with a concomittant discontinuity in the chemical potential ||. In 
sharp contrast, in NCCO, characteristic signatures of stripe order, e.g. the 1/8 anomaly 
and the NQR wipeout, appear greatly attenuated if not absent H, and the Fermi level 
shifts smoothly into the upper Hubbard band || . Here we show how the salient features of 
the most recent high resolution ARPES data in NCCO || can be understood in terms of 
the behavior of a uniformly doped Mott insulator. The doping dependence of the effective 
Hubbard U parameter adduced from the experimental spectra is in reasonable accord with 
various estimates, including a computation in which we screen the bare U via interband 
excitations. A quantum critical point (QCP) where the Mott gap closes is predicted just 
beyond optimal doping. Finally, we have carried out extensive first-principles simulations of 
photointensities in order to ascertain that the key spectral features discussed in this article 
(e.g. the appearance and growth of the (tc/2,tc/2) centered hole orbit with doping) are 
genuine effects of electron correlations beyond the conventional LDA framework and not 
related to the energy and k-dependencies of the ARPES matrix element || . 

Our analysis of the doping dependence of the Mott gap - perhaps better referred to 
as a pseudogap |||7J, invokes a mean field solution to the one band Hubbard model. 
The mean field approach provides a good description of not only the pseudogap in the 
ID charge density wave (CDW) systems flQfl , but also of the undoped insulator in the 



2D Hubbard problem, including the presence of the 'remnant Fermi surface' (rFS) (TT| 



seen in ARPES experiments [|12] and Monte Carlo simulations |13[]; the spin density wave 
(SDW) excitations are described by fluctuations about the mean field [14j]. The mean field 
theory correctly predicts a (stripe) phase instability (negative compressibility) associated 
with hole doping |Ll|; whereas assuming a finite second- neighbor hopping t' < 0, for electron 
doping the compressibility is positive [|nj, suggesting that a uniformly doped AFM state 
should be stable. The preceding considerations argue that the mean field would provide a 
reasonable model for discussing the doping dependence of the pseudogap. This viewpoint is 
further supported by mode-coupling calculations 0, which have been applied previously to 
approximately describe stripes in terms of a CDW in hole-doped cuprates fl7|Jl8| . 

Concerning methodology, we specifically consider the one-band Hubbard model where 
the Mott gap A is well known to arise (in the mean field SDW formulation Jl^j) from 
a finite expectation value of the magnetization at the wave vector Q = The 
self-consistent gap equation is 

i = uj2 — k —^ (i) 

k ^Ok 

where / is the Fermi function and 

Er = \^k + ^ ±E ok) (2) 

Here, the superscript c refers to the upper Hubbard band (UHB) and goes with the + sign 
on the right side while v goes with the - sign and the lower Hubbard band (LHB) = 
V ( e fc ~ e k+q) 2 + 4A 2 , and the independent particle dispersion is: = —2t(c x + c y ) —4t'c x c y , 
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Ci = cos kid. We study Eq. [T] self-consistently as a function of electron doping, assuming that 
the mean-field transition temperature corresponds to the experimentally observed pseudo- 
gap. U is thus treated as an effective parameter U e ff to fit the experimental data. The 
Green's function in the antiferromagnetic (AFM) ground state is given by 

k k 

where 

k 2 l 2E Q% h * 2 l 2E Q% ' 1 ; 

are the coherence factors. 

— * 

Using Eqs. 1-4, we discuss spectral density A(k, E) (given by the imaginary part of the 
Green's function) and the related energy dispersions for various doping levels. The FS maps 
are obtained by taking appropriate cuts through A(k, E). Although such theoretical maps 
do not account for the effects of the ARPES matrix element ||, these maps are relevant 
nevertheless in gaining a handle on the FS topology expected in the correlated system from 
ARPES experiments. In any event, we have also carried out first-principles simulations based 
on the conventional LDA picture where the ARPES matrix element is properly treated and 
the photoemission process is modeled including full crystal wavefunctions in the presence of 
the surface - see Refs. [0,0 f° r details of the methodology. 

Fig. p] displays the doping dependence of the spectral weight in the vicinity of the 
(tt/2, 7r/2) and (it, 0) points. Fig. |I|(a) shows that around (tt/2, it/2), with increasing doping, 
the spectral weight shifts rapidly towards the Ep as the gap between the LHB and UHB 
decreases. At x = 0.15, the gap is quite small and the LHB and UHB overlap. In con- 
trast, in the momentum region near (7r, 0), Fig. pjb), both LHB and UHB are more directly 
involved. Even at small doping levels (x = 0.04), Ep intersects the bottom of the UHB. 
With increasing electron concentration, the UHB moves to lower energies while the LHB 
shifts closer to the Ep as the gap decreases. These results are remarkably consistent with 
the corresponding ARPES data; Fig. 2a of Ref. shows a rapid movement of spectral 
weight at (tt/2, tt/2) from 1.3 eV to around 0.3 eV binding energy close to Ep. This effect, 
although somewhat less clear, is seen at (tt, 0) as well (Fig. 2b of Ref. PJ): the weight near 
Ep grows faster with doping, but the LHB is at a lower binding energy, and is less clearly 
resolved from the background. This doping dependence suggests that electrons first enter 
the UHB near (tt,0). 

The energy dispersions of Fig. 2 are useful not only in understanding the spectra of Fig. 
1, but also provide a handle on the FS maps expected in ARPES experiments. The Ep is 
seen in Fig. 2 to rise smoothly with respect to the UHB with increasing electron doping, 
whereas in the hole doped cuprates, the Ep gets pinned by the stripes near the mid-gap 
region over a large doping range. We should keep in mind that various bands do not possess 
the same spectral weight - this point is emphasized by depicting the coherence factors of 
Eq. 4 via the width of the bands in Fig. 2. As the Mott gap nearly collapses, the thick lines 
in Fig. 2(d) essentially present the appearance of the uncorrelated band structure. 

Fig. 3 illustrates the evolution of the FS corresponding to the band structure of Fig. 2. 
For low x, Ep lies near the bottom of the UHB and gives rise to electron pockets centered 
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at (tt, 0) and (0, tt) points in Fig. 0(a). The shadow segments of the FS are clearly evident in 
Fig. ^](b) where half of the electron pocket gains spectral weight at the expense of the other 
half due to the decreasing intensity in the magnetic Brillouin Zone (BZ) via the coherence 
factors; a weak imprint of the magnetic BZ boundary (the diagonal line connecting (0, tt) 
and (tt, 0) points) can also be seen. By x = 0.15, the gap shrinks considerably, and the 
Ep begins intersecting the LHB around (7r/2,7r/2). The FS now consists of three sheets: 
electron-like sheets near (tt, 0) and (0, it) and a hole-like sheet around (tt/2, tt/2)) separated 
by a residual gap located at the 'hot-spots' |^TJ along the BZ diagonal. Interestingly, in this 
doping regime, transport studies find evidence for two band conduction, and a change in the 
sign of the Hall coefficient [ 2"2"fl ; the hole pocket associated with the LHB can explain both 
these effects. 

We stress that the FS evolves following a very different route in the hole-doped case for 
the t — i! one band Hubbard model ]23| (neglecting stripes). At low doping levels, the FS 
consists of small hole pockets around [tt/2, tt/2) points. With increasing hole concentration, 
these pockets increase in size and merge to yield a large (tt, tt) centered FS satisfying the 
volume constraints of the Luttinger theorem. 

Figure 4 directly compares the theoretical FS maps against the corresponding experi- 
mental results. Here we have included a small second neighbor hopping parameter t" = O.lt 
in the computations in order to account for the slight shift of the center of the hole pocket 
away from (tt/2, tt/2) in the experimental data, even though the maps of Fig. 3 with only 
t — t' already provide a good overall description (after resolution broadening) of the measure- 
ments. The agreement between theory and experiment in Fig. 4 is remarkable: Both sets of 
maps show the electron pockets at low doping (x = 0.04); the beginning of the hole pocket 
at x = 0.10, together with the presence of shadow features around electron pockets and the 
appearance of intensity around the magnetic BZ; and finally, at x = 0.15, the evidence of a 
well formed hole pocket and a three-sheeted FS which has begun to resemble the shape of a 
large (tt,tt) centered hole sheet. Fig. 4 displays discrepancies as well (e.g. the hole pocket 
related feature appears double peaked in experiments, but not in theory, and other details 
of spectra around the Ep), but this is not surprising since we are invoking a rather simple 
one-band Hubbard model and the effects of the ARPES matrix element are missing in these 
calculations. 



The parameters used in the computations of Fig. 4 are as follows. |23] For the half-filled 
case, the values are identical to those used previously p5| : t = 0.326eV^, t'/t = —0.276, 
U = 6t. For finite doping x, the only change is that U is assumed to be x-dependent, 
U = U e ff(x), and t" = O.lt. The actual values of U e ff used (solid dots in Fig. 5(a)), are 
in reasonable accord with various estimates shown in Fig. 5(a) which are: the approximate 



value from Kanamori (arrow) for a nearly empty band , and Monte Carlo results p7,28 
for t' = (star). Finally, we have carried out a computation of screened U using U e ff = 
i7/(l+ < P > U), taking P as the charge susceptibility which includes only interband 
contributions with bare U = 6.75t. The U e ff so obtained for electrons (solid line) and holes 
(dashed) is shown P§| . 

Fig. 5(b) considers the behavior of the staggered magnetization, m^(x), for electrons 
and holes using the computed values of {7 e // for holes given by dashed line in Fig. 5(a). 
30| We see that m^(x) and hence the pseuodgap, which is proportional to m^(x), vanishes 



slightly above optimal doping, yielding a QCP. Notably, superconductivity near an AFM 
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QCP has been reported in a number of systems |3T|. However, the present case is different 
in that we have a mean-field QCP associated directly with short-range AFM fluctuations 
and the closing of the Mott gap. Since there is no interfering phase separation instability, 
the bulk Neel transition persists out to comparable, but clearly lower dopings, with T N — > 
near x = 0.13. 

We have carried out extensive first-principles simulations of the ARPES intensities in 
NCCO within the LDA framework for different photon energies, polarizations and surface 
terminations, in order to ascertain the extent to which the characteristics of measured FS 
maps could be confused with the effects of the ARPES matrix element || missing in the 
computations of Figs. 3 and 4. The computed FS maps including full crystal wavefunctions 
of the initial and final states in the presence of the Nd-CuCVNd-Cvterminated surface at 
16 eV are shown in Fig. 6 for two different polarizations, and are typical. The intensity 
is seen to undergo large variations as one goes around the (ir, ir) centered hole orbit, and 
to nearly vanish along certain high symmetry lines in some cases. Nevertheless, we do not 
find any situation which resembles the doping dependencies displayed in Figs. 3 and 4. It is 
clear that strong correlation effects beyond the conventional LDA-based picture are needed 
to describe the experimental ARPES spectra and that the t—t'—t" Hubbard model captures 
some of the essential underlying physics. 

In conclusion, this study indicates that the electron doped NCCO is an ideal test case for 
investigating how superconductivity arises near a QCP in doped Mott insulators, untroubled 
by complications of stripe phases. The doping dependence of Z7 e // adduced in this work 
has implications in understanding the behavior of the cuprates more generally since the 
pseudogap in both the electron and the hole doped systems must arise from the same Mott 
gap at sufficiently low doping. 
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FIGURES 



FIG. 1. Integrated spectral weight is shown near (a): (7r/2,7r/2), and (b): (w, 0) points, for 
four different doping levels, x = 0.0-0.15, as marked in (b). Domains in the Brillouin zone over 
which the spectral weight was integrated are shown in the insets. Fermi energy defines the energy 
zero in all cases. 



FIG. 2. Energy dispersions for various doping levels x. Energy zero defines the Fermi energy as 
in Fig. 1. Thickness of lines represents the spectral weights of various bands given by the coherence 
factors of Eq. 4. 



FIG. 3. Fermi surfaces corresponding to the band structures of Fig. 2. Maps are obtained by 
integrating the spectral density function (proportional to Im[G], Eq. 3) over an energy window of 
60 meV around Ep\ highs denoted by red and lows by blue. Experimental resolution effects are 
not included for the illustrative purpose of this figure. 



FIG. 4. Upper panels: Theoretical Fermi surface maps including resolution broadening for the 
t — t' — t" model (see text) for different doping levels x. Lower panels: Corresponding experimental 
maps after Ref. Q. Color scheme as in Fig. 3. 



FIG. 5. (a): U e ff (scaled by nearest neighbor hopping parameter t), and (b): staggered mag- 
netization mrf(x), vs doping x for electrons (solid lines) and holes (dashed lines). Filled circles give 
values used in the computations of Fig. 4 and are representative of the experimental ARPES data 
in NCCO. Values of U e ff given by Kanamori [gg] (arrow) and Monte-Carlo studies [§^J2^] (star) 
are shown in (a). 



FIG. 6. Theoretical FS maps in NCCO obtained via first principles simulations which include 
the effect of the ARPES matrix element but not of strong correlations for two different polarizations 
(given by the white arrows) at a photon energy of 16 eV. Color scheme as in Fig. 3. 
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